Radiative corrections to neutral pion-pair production 

N. Kaiser 

Physik-Department T39, Technische Universitat Miinchen, D-85747 Garching, Germany 



o 



Abstract 

We calculate the one-photon loop radiative corrections to the neutral pion-pair pho- 
toproduction process 7r~7 — > 7r _ 7r 7r . At leading order this reaction is governed by the 
chiral pion-pion interaction. Since the chiral tt + tt~ — > tt°tt° contact-vertex depends only 
on the final-state invariant-mass it factors out of all photon-loop diagrams. We give an- 
alytical expressions for the multiplicative correction factor R ~ a/2-rr arising from eight 
classes of contributing one-photon loop diagrams. An electromagnetic counterterm has 
£jy to be included in order to cancel the ultraviolet divergences generated by the photon- 

loops. Infrared finiteness of the virtual radiative corrections is achieved (in the standard 
! way) by including soft photon radiation below an energy cut-off A. The radiative correc- 

t^j- | tions to the total cross section vary between +2% and —2% for center-of-mass energies 

from threshold up to Tm^. The finite part of the electromagnetic counterterm gives an 
additional constant contribution of about 1%, however with a large uncertainty. 
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1 Introduction and summary 

The pions (tt + , 7r°, 7r~) are the Goldstone bosons of spontaneous chiral symmetry breaking in 
QCD. Their strong interaction dynamics at low energies can therefore be calculated system- 
atically (and accurately) with chiral perturbation theory in form of a loop expansion based 
on an effective chiral Lagrangian. The very accurate two-loop predictions [1] for the S-wave 
7T7r-scattering lengths, a = (0.220 ± 0.005)171^ and a 2 = (-0.044 ± O.OO^m" 1 , have been 
confirmed experimentally by analyzing the nir final-state interaction effects occurring in var- 
. ^ | ious (rare) charged kaon decay modes [2, 3, 4]. Electromagnetic processes with pions offer 
further possibilities to test chiral perturbation theory. For example, pion Compton scattering 
7r~7 — > 7T _ 7 allows one to extract the electric and magnetic polarizabilities (a n and (3 n ) of 
the charged pion. Chiral perturbation theory at two-loop order gives for the dominant pion 
polarizability difference the firm prediction a^ — (3 n = (5.7± 1.0) • 10 -4 fm 3 [5]. It is however in 
conflict with the existing experimental results from Serpukhov a^ — fa = (15. 6 ±7.8) • 10 _4 fm 3 
[6] and MAMI a n - fa = (11.6 ± 3.4) • 10" 4 fm 3 [7] which amount to values more than twice 
as large. Certainly, these existing experimental determinations of a, — fa raise doubts about 
their correctness since they violate the chiral low-energy theorem notably by a factor 2. It is 
worth to note that a recent dispersive analysis [8] of the Belle data for 77 — > 7t + tt~ gives the 
fit value a n — (3 n = 4.7 • 10 _4 fm 3 , compatible with chiral perturbation theory. 

In that contradictory situation it is promising that the ongoing COMPASS experiment 
[9] at CERN aims at remeasuring the pion polarizabilities, a n and (3 n , with high statistics 
using the Primakoff effect. The scattering of high-energy negative pions in the Coulomb-field 
of a heavy nucleus (of charge Z) gives access to cross sections for 7r~7 reactions through 
the equivalent photon method [10]. The consistent theoretical framework to extract the pion 
polarizabilities from the measured cross sections for (low-energy) pion Compton scattering 
7r _ 7 — > 7r~7 or the primary pion- nucleus bremsstrahlung process n~Z — > Ti^Z'y has been 
described (in one-loop approximation) in refs.fll, 12]. It has been stressed that at the same 
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order as the polarizability difference a w — j3 n there exists a further (partly compensating) pion- 
structure effect in form of a unique pion-loop correction (interpretable as photon scattering off 
the "pion-cloud around the pion"). In addition to these strong interaction effects, the QED 
radiative corrections to real and virtual pion Compton scattering 7r~7 ( *) — > 7T~7 have been 
calculated in refs.[12, 13]. The relative smallness of the pion-structure effects in low-energy 
pion Compton scattering [11] makes it necessary to include such higher order electromagnetic 
corrections. The COMPASS experiment is set up to detect simultaneously various (multi- 
particle) hadronic final-states which are produced in the Primakoff scattering of high-energy 
pions. The neutral pion production channel 7r~7 — > 7r~7r° serves as a test of the QCD chiral 
anomaly by measuring the 73-7T coupling constant = e/(47r 2 /^) = 9.72 GeV" 3 . For the 
two-body process 7r~7 — > 7r~7r° the one-loop [11, 14] and two-loop corrections [15] of chiral 
perturbation theory as well as QED radiative corrections [16] have already been worked out. 

The 7r~7 reaction with three charged pions in the final-state is used in the energy range 
above 1 GeV to study the spectroscopy of non-strange meson resonances [17] and to search for 
so-called exotic meson resonances [18]. The very high statistics of the COMPASS experiment 
allows it to continue the event rates with three pions in the final-state even downward to the 
threshold. The (differential) cross sections for the 7r~7 — > 3ir reactions in the low-energy region 
offer new possibilities to test the strong interaction dynamics of the pions as predicted by chiral 
perturbation theory. In a recent work [19] the production amplitudes for 7r _ 7 — > 7r~7r°7r° and 
7r~7 — > 7t + 7t~7t~ have been calculated analytically at one-loop order in chiral perturbation 
theory. It has been found that the next-to-leading order corrections from chiral loops and 
counterterms enhance sizeably (by a factor 1.5 — 1.8) the total cross section for neutral pion- 
pair production 7r~7 — > 7r~7r°7r°. By contrast the total cross section for charged pion-pair 
production 7T~7 — > 7t + 7t~tt~ remains almost unchanged in comparison to its tree-level result. 
This different behavior can be understood from the varying influence of the chiral corrections 
on the pion-pion final-state interaction (ix + it~ — >■ 7r°7r° versus 7r~7r~ — > 7r~7r~). 

The purpose of the present paper is to further improve the theoretical description of the 
7r~7 — > 37r reactions by considering the corresponding QED radiative corrections. We restrict 
ourselves here to the simpler case of neutral pion-pair production 7r~7 — > 7t~tt 7t , for which 
the number of contributing one-photon loop diagrams is limited to about a dozen. Another 
fortunate circumstance is that the (leading-order) chiral 7r + 7r~ — >■ it it contact-vertex factors 
out of all photon-loop diagrams and therefore the radiative corrections to 7r~7 — >■ 7r~7r°7r° can 
be represented simply by a multiplicative correction factor R ~ a/2n. Infrared finiteness of 
these virtual radiative corrections is achieved (in the standard way) by including soft photon 
radiation below an energy cut-off A. Taking A = 5 MeV, we find that the radiative corrections 
to the total cross section for 7T~7 — > 7r~7r°7r° vary between +2% and —2% for center-of-mass 
energies from threshold up to Im^. An electromagnetic counterterm (necessary in order to 
cancel all ultraviolet divergences generated by the photon-loops) gives an additional constant 
contribution of about 1%, however with a large uncertainty. The radiative corrections to the 
charged pion-pair production process 7r~7 — > 7t + 7t~tt~ can be roughly estimated to be a factor 
2 — 4 times larger, arguing that in this case twice as many charged pions are involved in virtual 
photon-loops and soft photon bremsstrahlung. 

2 Evaluation of one-photon loop diagrams 

In this section we calculate analytically the radiative corrections to the neutral pion-pair 
photoproduction process: 

7r-(pi)+7(fc,e) ^Tr-^ + Tr^gO + vr ^), (1) 
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as they arise from one-photon loop diagrams at order a. For a concise presentation of our an- 
alytical results it is convenient to introduce the following dimensionless Mandelstam variables: 

s = (pi + k) 2 /ml, t = (p 1 -p 2 ) 2 /ml , u= (p 2 -k) 2 /ml, (2) 

with m-K = 139.570 MeV the charged pion mass. In this (adapted) notation y/sm^ is the total 
center-of-mass energy of the process. We will also use frequently the linear combination: 

Y 1 = s + t + u-2 = {q l + q 2 ) 2 /ml , (3) 

related to the squared invariant mass of the produced neutral pion-pair. In the physical region 
the following inequalities hold: s > (1 + 2 v /7o) 2 , t < 0, u < 1 and Atq < X < (-y/s — l) 2 
where r = (m^o/m^) 2 = 0.93526 denotes the squared ratio between the neutral pion mass 
m w o = 134.977 MeV and the charged pion mass m n . 

Let us recall the dynamical description of the process 7r~7 — > 7r~7r°7r° at low energies 
[11, 19]. When choosing for the (transversal) real photon ^y(k, e) the Coulomb-gauge in the 
center-of-mass frame, the conditions e • p\ = e • k = imply that all diagrams for which the 
photon couples to the in-coming pion 7i~(pi) vanish identically. Furthermore, in the convenient 
parametrization of the special-unitary matrix-field U = \Jl — f 2 / f 2 + if ■ f/f n no ^An and 
2^An contact- vertices exist (at leading order). Under these assumptions one is left with one 
single -u-channel pole diagram in which the chiral 7r + 7r _ — > 7r°7r° contact-vertex is followed by 
a photon-pion coupling proportional to e • p 2 . 

The virtual radiative corrections to 7r~7 — > 7r~7r°7r° are obtained by dressing this tree 
diagram with a photon-loop in all possible ways (see Figs. 1-4). A fortunate circumstance 
is that the (leading order) chiral n + n^ — > 7r°7r° transition amplitude [(gi + q 2 ) 2 — m 2 { ,]/ f 2 
depends only on the 7r°7r° invariant mass and thus factors out of all photon-loop diagrams. 
The Coulomb-gauge (e-pi — e-k — 0) leaves the scalar product e-p 2 as the only possible coupling 
term for the external real photon. As a consequence of these features the radiative corrections 
due to photon-loops can represented simply by a multiplicative correction factor. Its real part 
which is only of relevance is denoted by R(s,t,u). We use dimensional regularization to treat 
both ultraviolet and infrared divergences (where the latter are caused by the masslessness of 
the photon). Divergent pieces of one-loop integrals show up in form of the composite constant: 

containing a simple pole at d — 4 and fi is an arbitrary mass scale. Ultraviolet (UV) and 
infrared (IR) divergences are distinguished by the feature of whether the condition for conver- 
gence of the ci-dimensional integral is d < 4 or d > 4. We discriminate them in the notation 
by putting appropriate subscripts, i.e. £uv and In order to simplify all calculations 

we employ the Feynman gauge, where the photon propagator is directly proportional to the 
Minkowski metric tensor g^ v . We can now enumerate the analytical expressions for R(s,t,u) 
as they emerge from the eight classes of contributing one-photon loop diagrams. 

The two diagrams of class (I) shown in Fig. 1 introduce the wavefunction renormalization 
factor Z 2 — 1 of the pion [12]: 

RP = -(tiR-tuv). (5) 



lr rhe inequality u < 1 follows immediately from the definition of u. In order to derive the sharper upper 
bound u < 0, one uses the relation for u written in eq.(21) and inserts y max = 1 and uj max — (s + 1 — Ar )/2y/s. 
In the end the condition r > 1/4 turns out to be crucial for u to take on only negative values. 
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Figure 1: One-photon loop diagrams (I) and (II) for neutral pion-pair production tt 7 — > 
7r~7r 7T°. Arrows indicate out-going pions. 
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Figure 2: One-photon loop diagrams (III), (IV) and (V). 



Diagram (II) involves the once-subtracted (off-shell) selfenergy of the pion and leads to the 
result: 

u + 1 , „ 



a 

Ti 



2u 



(6) 



Diagram (III) shown in Fig. 2 gives rise a constant vertex correction: 
while diagrams (IV) and (V) generate ^-dependent vertex corrections: 



(7) 
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It is astonishing that the last four contributions R^ + R( m > + R^ + i?( v ) = sum to zero. 

The reducible u-channel pole diagram (VI) shown in Fig. 3 includes a photonic vertex 
correction around the 2n° emission vertex. One finds for its contribution to the (real) R- 
factor the following result: 



R( yi) = «. 

2tt 



iw + 1 + ln(l - u) + 4/ - 



4, v/S^I+v^ 
In — 



+ s + t + 



■in ~ v ' \ >; 2 

u-7\ f 1 , hilar 1 + £(z - 1)| - ln(l - u) 
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Jo 
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(10) 
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The integrand of the principal- value integral dx has simple poles at x± = [H + 1 — u ± 
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Figure 3: One-photon loop diagrams (VI) and (VII). 



inside the unit-interval < x < 1. Due to this property an accurate numerical treatment of 
the principal- value integral in eq.(10) (in combination with further three-body phase space 
integrations, see eq.(20)) is easily manageable. We have also checked the Feynman parameter 
representation of the last loop integral in eq.(10) against its dispersion relation representation: 



1 In {x" 1 + E(x - 1)| - ln(l - u) _ r°° dx 

1 + (u-1)j; + Si(i-1) h x-S J( x + l-u) 2 -4x 



x In 



y/x(x - u - 3) + JJx - 4) \(x + 1 - u) 2 - Ax] 



2(1 -it) 



(11) 



where the imaginary part on the right hand side has been calculated via the Cutkosky cutting 
rule. The analytical continuation of the last logarithmic term in the first line of eq.(10) for 

'E/2). 



< E < 4 is J 4/E - 1 arcsin( 

In the sum + i?( VI ) obtained so far the ultraviolet divergences £uv do not cancel and 
the remaining classes of diagrams (VII) and (VIII) are actually ultraviolet convergent. In a 
non-renormalizable effective field theory like chiral perturbation theory such a behavior of the 
radiative corrections is generic. In order to eliminate all ultraviolet divergences from photon- 
loops additional electromagnetic counterterms have to be introduced [20]. The black square 
in the right tree diagram of Fig. 4 symbolizes the pertinent electromagnetic counterterm for 
7r + 7r~ — > tt°tc° scattering. It gives rise to the following constant contribution to the .R-factor: 



R (ct) 



3a 
2^ 



{£uv + k) 



(12) 



where k denotes the finite part of the electromagnetic counterterm which remains after cancel- 
ing the ultraviolet divergence £uv- A numerical estimate of k will be given in section 4. For the 
sake of completeness we quote the expression for the on-shell 7r + (g + )+7r~(g_) —> 7t°('?i)+71" (q , 2) 
scattering amplitude with inclusion of radiative corrections [20]: 
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(13) 



where E = (g x + q 2 ) 2 /m 2 r and f n = 92.4 MeV denotes the pion decay constant. In the order 
given the terms in the square bracket correspond to the pion wavefunction renormalization 
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Figure 4: One-photon loop diagrams (VIII). The black square in the right tree diagram (ct) 
symbolizes the electromagnetic counterterm for tt + 7t~ — > 7r°7r° scattering. 



factor Z 2 — 1, the electromagnetic counterterm, and the one-photon exchange contribution. 
Note that we have used the (concise) spectral function representation for the infrared divergent 
(scalar) loop integral involving one photon and two pion propagators. If an infinitesimal photon 
mass m 7 is introduced as an (alternative) infrared regulator the infrared divergence £m is to 
be identified with the logarithm In (m^/m^). 

Next, we come to the irreducible s-channel diagrams of class (VII) shown in Fig. 3. Since 
the contribution from the left diagram (involving the two-photon contact-vertex) gets (partly) 
canceled by a term from the right diagram it is advantageous to specify only their total 
contribution to the i?-factor. After reducing the loop integrals with four propagators one 
obtains the following result: 
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(14) 



Note that the last denominator has no poles in the physical region, since s — 1 + S(x — 1) > 
s — 1 — E > s — 1 — (y/s — l) 2 = 2(y / i — 1) > 0. By taking the absolute magnitude of 
the arguments of logarithms one gets directly a suitable representation of the only relevant 
real part. It is a fortunate circumstance that the Feynman-parameter representation of loop 
functions leads to expressions which can be handled easily numerically in the physical region. 

Finally, we come to the irreducible w-channel diagrams of class (VIII) shown in Fig. 4. The 
contribution from the left diagram gets completely absorbed by a term from the right diagram. 
The resulting contribution to the (real) i?-factor includes an infrared divergent term with a 
non-trivial t-dependence and after putting all pieces together it reads: 
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(15) 



fi 



with the abbreviations 



w = \{^-{^i)i h ± = \{t±^ r ^t). (16) 
One observes that the term proportional to D(t) — -D(S) drops out in the sum 

i? (VII) +jR (VIII) 

and therefore we do not need to specify it. Li 2 (u>) = X^li n~ 2 w n = w dx[x(x — w)]^ 1 \nx 
denotes the conventional dilogarithmic function. Several of the results derived in section 3 of 
ref. [13] have been useful in order to obtain the expression for i?( vm ) written in eq.(15). 



3 Infrared finiteness 



In the next step we have to consider the infrared divergent terms proportional to present in 
eqs.(5,15). At the level of the measurable cross section these get eliminated by contributions 
from (undetected) soft photon bremsstrahlung. In its final effect, the (single) soft photon 
radiation off the in- or out-going n~ multiplies the tree-level differential cross section for 
7T~7 — >■ 7r~7r°7r° by a (universal) factor [12, 13]: 
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(17) 



which depends on a small energy cut-off A. Working out this momentum space integral by the 
method of dimensional regularization (with d > 4) one finds the following contribution from 
soft photon emission to the i?-factor: 
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(18) 



with the abbreviation W = (s + 1 — Ei) 2 — 4s[l — tx(l — x)]. In order to simplify the last 
term in eq.(18) we have made use of the relation E = s + 1 — 2ajy/s, where oj denotes the 
center-of-mass energy of the out-going negative pion 7r~(p 2 ) divided by m n . Note that the 
terms beyond those proportional to ln(m„-/2A) — (,ir are specific for the evaluation of the soft 
photon correction factor i? so ft in the center-of-mass frame with A an infrared cut-off therein. 

In order to present a concrete example we have evaluated the complete radiative correction 
factor R at the threshold in the isospin limit: s t h = 9, t t h = ~ 4/3, w t h = —5/3, S t h = 4, 
Uth = 1. In this case one gets numerically: 
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(19) 



where the terms in the curly bracket correspond in the order written to virtual photon-loops, 
the electromagnetic counterterm, the universal soft photon contribution, and the soft photon 
contribution specific for imposing an infrared cut-off via \l\ < A in the center-of-mass frame. 
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Figure 5: Radiative corrections to the total cross section for neutral pion-pair photoproduction 
7T~7 — > 7r~7r°7r° as a function of the center-of-mass energy ^/sm n . The infrared cut-off for soft 
photon emission has been set to the value A = 5 MeV. 

4 Results: radiative corrections to cross sections 

After inclusion of radiative corrections the total cross section for neutral pion-pair production 
7r~7 — > 7r~7r°7r° depends also on the infrared cut-off A for undetected soft photons. We multiply 
the squared tree-level amplitude by 1 + 2R(s, t, u, A) and integrate over the three-pion phase 
space. Applying the usual flux and symmetry factors the total cross section reads: 
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(20) 



with o; max = (s + 1 — Ar )/2y/s the endpoint energy of the out-going n divided by m n . Using 
the relations S = s + t + -u — 2 = s + l— 2ujy/s and: 



l + ^^(w-yVw 2 -l), t = 2-(s + 1)4= + ^^yVw 2 - 1 , 



(21) 



valid in the center-of-mass frame the whole integrand in eq.(20) becomes a function of u and 
the directional cosine y. 

Fig. 5 shows in percent the radiative corrections to the total cross section for neutral pion- 
pair production 7r~7 — > 7r~7r°7r° as a function of the center-of-mass energy y/sm n . The dashed- 
dotted and dashed curves display the separate contributions from soft photon bremsstrahlung 
and virtual photon-loops. In each case the radiative correction is calculated as the ratio of 
the shift in cr to t(s, A) induces by the respective i?-factor divided by the tree-level cross section. 
As in ref . [12] the infrared cut-off A for undetected soft photons has been set to A = 5 MeV, 
a value which seems appropriate for the COMPASS experiment. The full line in Fig. 5 shows 
the complete radiative corrections. One observes an almost linear decrease which ranges from 
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Figure 6: Radiative corrections to the 7r°7r° mass spectra for neutral pion-pair production 
7r~7 — > 7r~7r°7r° as a function of the 7r°7r° invariant mass v^S^tt- The numbers on the curves 
correspond to y/s. 

+ 1.6% at threshold to —1.6% at the center-of-mass energy y/sm^ = lm^. An interesting 
feature is that the positive radiative corrections from the virtual photon-loops get gradually 
reduced and turned into negative values by the soft photon contributions. 

The finite part of the electromagnetic counterterm k shifts the radiative corrections (dis- 
played by the full curve in Fig. 5) by a constant amount of 3ak/n = 0.7% • k. In order to give 
an estimate for k we exploit the elaborate result of ref. [21] for the pionium decay amplitude 
(a — + £■ Guided by eq.(13) we identify (a/27r)(3A; + 1) with the ratio e chn / (a — 02), 

where e elm is the electromagnetic correction to the pionium decay amplitude. Subtracting from 
e = (6.1±1.6)-1CT 3 the contribution 4.8-10~ 3 due to the (charged and neutral) pion mass differ- 
ence (see eqs. (4.28,4.29) in ref. [21]) one gets e chn = (1.3±1.6)-10 -3 . Together with the leading 
order expression for the 7T7r-scattering length difference a — a 2 = 9m n /(327rf^) = 0.204m" 1 
one arrives at the estimate k = 1.5 ±2.2 for the electromagnetic counterterm. Its central value 
implies a constant shift of the radiative corrections to 7r~7 — > 7r~7r°7r° by about 1.0%. The 
large errorbar of k = 1.5 ± 2.2 introduces at the same time a wide errorband to the full curve 
in Fig. 5. Still an allowed option is to neglect to electromagnetic counterterm, setting k — 0. 

Finally, we consider radiative corrections to more exclusive observables. An obvious can- 
didate is the 7r°7r° mass spectrum da/dmoo with m o = \/2m ff the 7r°7r° invariant mass. The 
differential cross section do-/dm 00 is obtained by omitting the (iw-integration in eq.(20) and 

applying the normalization factor m~ 1 \J'E/s. Fig. 6 shows in percent the radiative corrections 
to the 7r°7r° mass spectrum for neutral pion-pair production 7r~7 — > 7r~7r°7r° as a function of 
the 7r°7r° invariant mass V^m,. The numbers (4,5,6,7) on the four rising curves correspond 
to y/s, the total center-of-mass energy divided by m^. The electromagnetic counterterm k 
shifts again the whole pattern by a constant 3ak/n. 

In summary, we find that the radiative corrections to neutral pion-pair production 7r~7 — >■ 
7r~7r°7r° are comparable in size to those for pion Compton scattering 7r~7 — > 7r~7 [12]. 
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